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The thermal boundary resistance of Si/Ge interfaces as been determined using approach-to- 
equilibrium molecular dynamics simulations. Assuming a reciprocal linear dependence of the ther¬ 
mal boundary resistance, a length-independent bulk thermal boundary resistance could be extracted 
from the calculation resulting in a value of 3.76x 10“® m^K/W for a sharp Si/Ge interface and ther¬ 
mal transport from Si to Ge. Introducing an interface with finite thickness of 0.5 nm consisting of 
a SiGe alloy, the bulk thermal resistance slightly decreases compared to the sharp Si/Ge interface. 
Further growth of the boundary leads to an increase in the bulk thermal boundary resistance. When 
the heat flow is inverted (Ge to Si), the thermal boundary resistance is found to be higher. From 
the differences in the thermal boundary resistance for different heat flow direction, the rectifica¬ 
tion factor of the Si/Ge has been determined and is found to significantly decrease when the sharp 
interface is moderated by introduction of a SiGe alloy in the boundary layer. 


I. INTRODUCTION 

Thermoelectricity as alternative energy production has 
gained increased interest in recent years and substantial 
research has emerged to increase the hgure of merit ZT 
of thermoelectric materials.^ilf^ The figure of merit de¬ 
scribes the efficiency of a thermoelectric material and 
depends on its Seebeck coefficient and its electric and 
thermal conductivity. A common approach to achieve a 
high ZT is by minimization of the thermal conductivity 
of semiconducting ma terials and preserving good detailed 

conduction properties 

Introduction of impurities to pristine bulk materi¬ 
als by alloying has been shown to be a promising way 
to minimize the material thermal conductivity n. In 
fact, the thermal conductivity of Sia;Gei_a; alloys is re¬ 
duced remarkably compared to their bulk counterparts 
already at very small impurity concentrations x.^SHoigini- 
ilar e ffect s have been shown in Bi 2 Te 3 , Bia;Sb 2 _a:Te 3 and 
PbTe.l^ Scattering of short-wavelength phonons on im¬ 
purity atoms is responsible for the decrease of n while 
thermal transport by mid- and long-wavelength phonons 
remains unperturbed by atomistic defects.^ It is therefore 
essential to suppress the propagation of phonons with 
longer wavelength for achieving a further reduction of k. 
A possible way to realize this is by nanostructuring of 
the material confirmed by several experimental and the¬ 
oretical studies. 

Recently, introduction of nanoscopic holes in a SiGe 
bulk alloy simulated by molecular dynamics has been 
shown to drastically decrease the number of low- 
frequency (long-wavelength) phonons resulting in a re¬ 
markable reduction of the thermal conductivity in the 
SiGe alloy.^ Furthermore, a lower thermal conductiv¬ 
ity in Si/Ge-based materials was observed in either su¬ 
perlattices or nanowires and nanodots depending on the 
dimension of the nanostructures.^ Similar results have 
been shown in experimental studies where the thermal 
conductivity of bulk Si and bulk SiGe alloys was reduced 


significantly when nano crystall ine structures have been 
generated by ball milling.EESIill The decrease of k in such 
nanostructured materials results from increased phonon 
scattering at the interfaces introduced by nanograins, 
nanowires or superperiodicity. 

The effect of interfaces to thermal conductivity is de¬ 
noted as thermal boundary resistance (TBR), also known 
as Kapitza resistance.li^ Yet, a complete understand¬ 
ing of the scattering properties of the interface affect¬ 
ing the TBR and an accurate prediction of such is still 
a matter of investigation. Several models have been 
proposed for the determination of TBR such as the 
diffuse and the aco ustic mismatch model (DMM and 
AMM, respectively) However, within these models 
assumptions of the scattering processes are made a pri¬ 
ori. Phonon interface scattering, for example, is assumed 
to be elastic in both methods. Furthermore, the DMM 
considers only diffuse scattering while in the AMM dif¬ 
fuse scattering is neglected and only specular scattering 
is assumed. The TBR values obtained from these models 
are in general lower than in experiments; for instance, 
in the example of Si/Si 02 interfaces experimental values 
range from 2xI0“® to 2xI0“® m^K/^^SlH 20 ].^J^J^g theory 
predicts a TBR of 2.4-3.5 xI0“® m^K/W.^ 

Another possible way to predict TBR is based on 
molecular dynamics (MD) simulations. Determination of 
the TBR of Si/Ge interfaces using non-equilibrium MD 
simulations resulted in a value of 1.26 x I0“® m^K/W at 
a process temperature of 300 K.I^In this study, the ther¬ 
mal boundary resistance has been determined at a finite 
simulation cell size of 20 nm which can significantly af¬ 
fect the thermal transport properties as has been shown 
for the length dependent thermal conductivity of bulk 
materials.®!^ A higher value (2.72-3.17 xIO"® m^K/W) 
has been reported by Landry et al. which was reason¬ 
ably close to the TBR predicted by DMM (2.4 xI0“® 
m^K/W).^ Moreover, non-equilibrium MD simulations 
have recently been used to determine the TBR in SiGe 
nanowires.^ 
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The TBR can be affected by the direction of the heat 
flux as a result of the asymmetry in thermal conduc¬ 
tivity of the bulk materials Si and Ge which is dehned 
in the rectification factorP^ Depending on the mass and 
lattice mismatch, the rectification of a two segment bar 
can change signihcantlyP^ Therefore, the quantification 
of rectihcation is also of great interest in the characteri¬ 
zation of the TBR of a certain material pair. 

Here, we have used approach-to-equilibrium MD 
(AEMD) simulations to estimate the TBR of Si/Ge inter¬ 
faces. First, we have derived an expression for the length- 
independent TBR at the theoretical limit of an infinite 
sample length (section |IIB ). Based on this expression. 


the TBR of a sharp Si/Ge interface has been evaluated 
for heat flux both from Si to Ge and vice versa (section 
Furthermore, the effect of boundary thickness has 


been in vestiga ted for a boundary thickness up to 2 nm 
(section IIIB). Finally, the rectification of the interface 
has been determined as a function of the boundary thick¬ 
ness by comparison of the TB R for heat flux from Si to 
Ge and from Ge to Si (section IIIG). 


II. METHODS 


A. AEMD simulations 


Molecular dynamics simul ations have been performed 
using the lammps code.l^^l^Interatomic forces have been 
described applying the Tersoff pair potentiaP^^ which has 
been shown to represent reasonably well the mecha nical 
and thermal properties of Si and Ge materials 

The overall thermal conductivity of the simulated sys¬ 
tems has been determined using the AEMD method. 
Originally, the AEMD method was developed under the 
assumption that k is uniform in the direction of the heat 
transport. This is obviously not the case for the here 
investigated systems. However, from considerations of 
energy conservation, we can argue that heat transport 
between two homogeneous systems, where k is stepwise 
uniform, can indeed be described with an overall thermal 
conductance coefficient according to 


q = hc‘ at 

, . La (1) 

with he = — + R -\ - 

Ki K2 

Assuming such an overall, uniform thermal conduc¬ 
tance coefficient for the differential heat transport, the 
same procedure and equations of the AEMD approach 
for homogeneous systems can be applied for heat trans¬ 
port through interfaces with stepwise uniform thermal 
conductivities, supporting the use of the latter also for 
heterogeneous systems where an overall constant thermal 
conductivity coefficient he can be applied. This assump¬ 
tion has been used previously for the calculation of the 
thermal conductivity and boundary resistance of Si/Si 02 


interfaces.!^ Furthermore, we anticipate a result exten¬ 
sively described in section HI A[ namely: for a sharp in¬ 
terface, AEMD results are in excellent agreement with 
other theoretical predictions. 

Using the AEMD method, the simulation cells are 
firstly divided into two regions with comparable length in 
the direction of thermal transport. One of these two com¬ 
partments is equilibrated at high temperature (Th=400 
K), the other compartment at low temperature (Tc=200 
K) using velocity rescaling. This creates an initial step¬ 
like function of the temperature along the sample length 
in z-direction.l^ Next, the evolution of the average tem¬ 
perature in the hot (Th) and cold (Tc) reservoir has been 
recorded during a transient regime towards equilibrium of 
microcanonical evolution. Based on Fourier’s theorem of 
thermal transport and the given step-like initial temper¬ 
ature profile, the evolution of the temperature gradient 
(AT = Th - Tc) follows 


OO 

AT(t) = ;^C„e-“"«*, (2) 

n—1 

where R = — is the thermal diffusivity with the density 
p of the material and its specific heat Cy. This expression 
is fitted to the temperature gradient obtained from the 
simulations to determine k. More det ails on the method¬ 
ology can be found elsewhere. d*^ l ^^ l ^^ l 

Creation of Si/Ge interfaces A schematic image of 
a Si/Ge simulation cell is represented in Fig. The 
Si/Ge crystals are oriented with the crystallographic 
(001) plane orthogonal to the heat flux (in-plane direc¬ 
tion). To account for periodic boundary conditions in 
the in-plane direction, the in-plane lattice spacing has to 
be equal for both Ge and Si sections. Motivated by ex¬ 
perimental studies,!^ where SiGe alloys have been grown 
on pure Ge, we have imitated growth of Si on crystalline 
Ge where the equilibrium lattice parameter of Ge (ace.o) 
has been adopted for the in-plane lattice spacing of Si. 
The out-of-plane lattice parameter of Si (asi,_L) has been 
determined from the elastic properties of Si according to 


ai’ = 


1 - 2 


(Ch2 

\Cu 


Si 



( 3 ) 


with Cii and C 12 of 142.54 and 75.38 GPa, respec¬ 
tively, resulting in asi,_L=5.1913 A. The elastic constants 
Cii and C 12 applying for the Tersoff potential have been 
determined from the second derivatives of the total en¬ 
ergy with respect to deformation. Pseudomorphic Si 
samples with these properties will be referred to as p- 
Si in the following. 

In samples with a finite interface thickness dj, several 
atomic layers of a Sia;Gei_a; alloy have been added at 
any interface. The Si concentration x in the finite inter¬ 
faces has been gradually increased every two atomic lay¬ 
ers from the pristine Ge section to the pure p-Si section. 
The interfacial spacing has been determined equivalent 
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FIG. 1. Schematic representation of a Ge/Si interface (in 
the (001) crystallographic plane) with boundary thickness dj. 
The lattice spacing in in-plane direction (orthogonal to ther¬ 
mal transport) of Si has been adopted from the equilibrinm 
lattice spacing of Ge (aGe,o=5.6567 A). The out-of-plane lat¬ 
tice spacing of Si has been calculated from the elastic con¬ 
stants (asi,x=5.1913 A). 


to p-Si based on eq. ^ and has been adapted for each 
Si concentration x. An interfacial spacing of 0.5, 1 and 2 
nm has been simulated, corresponding to a total of four, 
eight and 16 atomic layers, respectively. 



1/L^ [nm'^] 

FIG. 2. Inverse of the thermal conductivity 1/k as a function 
of the inverse sample length 1/Lz of crystalline Ge (c-Ge, dia¬ 
monds), crystalline Si (c-Si, triangles up) and pseudomorphic 
Si (p-Si, triangles down). The bulk thermal conductivity Kcx, 
of these materials has been approximated based on a linear 
relationship between 1/k and l/Lz- It resulted in 93.3, 233.4 
and 204.3 W/mK for c-Ge, c-Si and p-Si, respectively. 


B. Determination of thermal boundary resistance 


The overall thermal conductivity Kail of a heteroge¬ 
neous system, such as the heterostructure shown in Fig. 
^ respectively its thermal resistance (i?aii = ) j can be 

aescribed as a connection of series of resistances. For the 
systems calculated here (Fig. this can be expressed 
according to 


i?aii = ^ = + 2Ri (4) 

Kail KGe(RGe) Kp-Si(Tp-Si) 


where Ace and Asi are the total simulation cell 
length and the length of the crystalline Ge and the pseu¬ 
domorphic Si part, respectively (see Fig. [^. Kce and 
Kp-Si represent the thermal conductivities of crystalline 
Ge and p-Si at their cell length Lqb and Lgi, respec¬ 
tively. The thermal boundary resistance Rj enters in the 
equation twice due to the periodic boundary conditions 
of the simulation cell. 

Thermal conductivity of k (L) of Ge and Si Following 
eq. the length-dependent thermal conductivities of 
the pure systems Kce and Kp.si have to be known for the 
determination of i?aii. This has been done using t he stan- 
dard AEMD approach for a homogeneous systerEpSUZIll] 
as described above. Therefore, the sample length of these 
systems has been varied from 100 nm (~200ao ) to pm 
( 2 OOO 00 ). The behavior of 1/k to 1/Az has been approx¬ 
imated by a linear function ^ ^ ^1 + > where Koo 

is the bulk thermal conductivity and A can be defined 
as characteristic length of the phonon transport. The 
linear approximation is a common way to describe the 
dependency of k on the sample length and is approved 
for systems where phonon properties are approximated 


TABLE I. Optimized parameters Koo and A, 
length-dependent thermal conductivity (^ = 
of crystalline Ge, Si and pseudomorphic Si. 


describing the 

±(i + A)) 




LmKJ 


A [nm] 


c-Gi 93.3 ± 3.5 430 ± 24.4 
c-Si 233.4 ± 13.4 742 ± 72.4 
p-Si 204.3 ± 20.6 855 ± 108 


well by an average value.l^SHH It has been shown to give 
reasonable estim ations o f the bulk thermal conductivity 
in Si/Ge systems .I22l2lll2l 

The length-dependent thermal conductivity of crys¬ 
talline Ge, Si and pseudomorphic Si is shown in Fig. 
with bulk thermal conductivities estimated to be 93.3, 
233.4 and 204.3 W/mK, respectively (Table]!]). The strain 
applied on the pseudomorphic Si as a result of the non¬ 
equilibrium in-plane lattice spacing thus reduces the ther¬ 
mal conductivity of the crystalline Si by 12.4 %. 

Overall thermal conductivity The overall thermal 
conductivity is a summation of several phonon transport 
effects (eq. (©)• Accordingly, the assumption of average 
phonon properties as in the case of the pure Ge and Si 
materials is unfounded. To account for the non-linear 
effects in Kaii, the Taylor expansion of 1/k to 1/Az has 
been extended to the second order term. 

— = —+ ^ + ( 5 ) 

Kail Kall^cso \ J 

In fact, a linear approximation results in a negative 
bulk thermal conductivity Kaii,oo, whereas the fitted sec- 
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FIG. 3. Total thermal conductivity K^ii as a function of the 
sample size Lz of Si/Ge interfaces with a sharp interface. 
Thermal transport has been simulated from Si to Ge (dia¬ 
monds) and vice versa (circles). The behavior of k to Lz 
has been approximated with a 2“'^ order polynomial function 
(dashed line). 


ond order function approximates the calculated values of 
Kali perfectly as shown in Fig. |^for sharp interfaces. 

TBR estimation In order to estimate the thermal 
boundary resistance the effects of a finite simulation cell 
have to be eliminated. Therefore, eq. Q is rewritten as 


Kali KQe {^Ge) Kp_Si {Lp-Si) 


where the length of the Ge and Si part is expressed 
with respect to the total cell length as ace = and 

asi = respectively. 

At infinite sample size (L^, —^ 00 ), the length- 
dependent Kali, KGe and Kp_si Converge to their bulk 
thermal conductivities Kaii.oo KGe,oo and Kp_si,oo, respec¬ 
tively, and the last term ^ of eq. (|^ vanishes. With 
these assumptions, the bulk thermal conductivities are 
related according to eq. 0. 


1 _ QlGe QlSi 

^all,oo ^Ge,oo ^p-Si,oo 


( 7 ) 


With the linear approximation of I/k( 1/L) in c-Ge 
and p-Si and the quadratic behavior of 1/Kaii (1/T) for 
the overall thermal conductivity (eq. (©), @ can be 
rewritten as follows 


Rj — - Tci X Lz + C2 + C3 X — 


where Ci = 


C2 = 


I 

^all,oo 

A. 


^all 


^all 


,00 


O^Ge _ QSi 

^Ge,oo ^p-Si,oc 

Agb _ Ap_si 

^Ge,oo ^p-Si,oc 

Mall 


G. = 


^all,oo 


( 8 ) 


FIG. 4. Thermal boundary resistance Ri as a function of the 
sample size Lz- The behavior of Rj to Lz follows a linear 

reciprocal trend Ri — Ri ,00 


As a result, the thermal boundary resistance this sim¬ 
plifies to 


i;, = i(G + C3xT)=s,„(i + ^) ,9) 

where i?/,oo is the converged TBR at bulk conditions 
{Lz 00 ) and A/ is the characteristic phonon length of 
the interface. To determine the parameters of Rj {Lz), 
Ri has been calculated at each data point from the over¬ 
all thermal conductivity Kaii (eq. ©) and previously 
determined length-dependent k of the pristine Ge and p- 
Si materials. The obtained values have then been fitted 
to eq. 


III. RESULTS AND DISCUSSION 
A. TBR at sharp interface 

In Fig. 1^ the calculated AEMD values of Ri as a 
function of the simulation cell length Lz together with 
their fits by eq. are shown for Si/Ge interfaces with 
various boundary thickness dj where thermal transport 
has been simulated from Si to Ge. 

The bulk TBR Rj^oo of a sharp interface has been 
calculated to be 3.76x10“® m^K/W (Table [ill)■ Several 
theoretical methods have been used in a previous study 
for the estimation of the TBR.I^ Among others, non¬ 
equilibrium molecular dynamics simulations have been 
used where the thermal resistance is calculated from the 
temperature gradient across the interface. Using this 
approach, the thermal boundary resistance at 300 K 
resulted in 3.1x10“® m®K/W. Despite the differences 
in the applied methods in ours and the previous work 
(which include a different lattice spacing, a different pair 
potential and a different approach to calculate Rj), this 
is in very good agreement with our results. Furthermore, 
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the value obtained from both MD-based methods are 
in reasonable agreement with the TBR calculated from 
the theoretical, Laundauer-like expressions for phonon 
scattering (3.0x10“® m®K/W) and the results using the 
DMM (2.4x10“® m®K/W(P^ where approximations of 
the nature of phonon scattering (such as elastic scatter¬ 
ing) have to be applied. The AEMD method used here 
thus provides a robust approach to reliably estimate the 
TBR without the necessity of approximations for phonon 
scattering and verifies the applicability of a mean ther¬ 
mal conductivity which is assumed to be uniform in the 
direction of the heat transport. 

When the heat flow is inverted, i.e. the Ge region is 
initially equilibrated at T/i=400 K and p-Si at Tc=200 K, 
the thermal boundary resistance increases to 5.78x10“® 
m®K/W (Table [n|. Similar results have been shown pre¬ 
viously for the heat transport in SiGe nanowires.^^ Using 
non-equilibrium MD simulations, the heat transport has 
been shown to be higher for a heat flow direction from Si 
to Ge than vice versa. 

The increase by ~ 50 % can be explained by the dif¬ 
ferent thermal conductivities of the pristine Ge and p-Si 
materials. Silicon has a higher thermal conductivity than 
Ge. Accordingly, the inflow of heat from the p-Si region 
to the interface is facilitated over its outflow through the 
Ge region which leads to a higher heat pressure at the 
Si/Ge interface. A higher heat pressure in this sense is 
comparable to a higher temperature of the system. This 
is comparable to previous results where the TBR has 
been shown to decrease with increasing process temper¬ 
ature from 300 to 1000 K.I^Sl 


B. Effect of interface thickness on TBR 

The morphology of the Si/Ge interface has been var¬ 
ied by changing its thickness dj in the range 0 nm< dj < 
2 nm. At a finite thickness, the interface has been con¬ 
structed of a Sia;Gei_a; alloy in which the Si concentration 
X has been increased gradually from the Ge region to the 
Si region as described in section [Ilj 

Switching from a vanishingly thin interface (di = 0 
nm) to a finite thickness of 0.5 nm results in a slight 
decrease (Ai?/ = 0.1 x 10“® m®K/W, see Fig. [^and Table 
0 of the TBR in the case of thermal transport from p-Si 
to Ge. This effect is more pronounced when the heat flow 
is inverted (Ai?/ = 1.1 x 10“® m®K/W). The introduction 
of an interface with finite thickness reduces the lattice 
and mass mismatch between the pure materials and the 
SiGe alloy. As a result, the TBR is reduced. 

Further boundary growth, however, results in a steady 
increase of the thermal boundary resistance. For thermal 
transport from p-Si to Ge, it raises up to 7.18 xl0“® 
m®K/W at a boundary thickness of 2 nm (Table 0- 
This is in agreement with previous results of the TBR in 
SiGe nanowires calculated by non-equilibrium MD sim¬ 
ulations showing TBR to increase when the Si/Ge inter¬ 
faces thickness increases from abrupt to 5 and 15 nm.l^ 



dj [nm] 

FIG. 5. Bulk thermal boundary resistance R/,oo as a function 
of the interfacial thickness dj for thermal transport from Si 
to Ge (diamonds) and from Ge to Si (circles). 

When the heat flow direction is inverted, it reaches 8.37 
xl0“® m®K/W. The increase in TBR with increasing 
thickness can be explained by increased phonon scatter¬ 
ing inside the boundary layer. For a sufficiently large 
boundary, this effect dominates the reduced lattice and 
mass mismatch and the TBR exceeds the value of the 
sharp interface. In the case of thermal transport from Si 
to Ge this is already the case at di=l nm while for heat 
flow from Ge to Si, it is observed at di =2 nm. 

The characteristic phonon length A/ (eq. |^, on the 
other hand, is hardly affected from the heat flow direc¬ 
tion. Except for the sharp interface. A/ (d) is within 
10% for both directions of the heat flow (Table 0. It 
decreases monotonously from 78.1 (76.3) to 26.2^7.8) 
nm with increasing boundary thickness from 0.5 to 2 
nm for a heat flow from p-Si to Ge (Ge to p-Si). The 
characteristic length does not directly correspond to any 
phonon property, however, it gives an indication. Here, 
its decrease with increasing boundary thickness can be 
regarded as a decrease in the phonon mean free path, 
presumably resulting from enhanced phonon scattering 
inside the boundary layer. 

C. Rectification 

Thermal rectification can be understood in analogy to 
an electrical diode. A material is defined as a thermal rec¬ 
tifier if the magnitude of the heat flow is different when 
inverting the heat flow direction.!^ Several expressions 
have been used recently to specify the thermal rectifi¬ 
cation. Most commonly it is defined as the difference 
of the magnitude of back ward and f orward heat flow di¬ 
vided by the smaller Within the AEMD 

framework, the magnitude of the heat flow is not directly 
calculated. However, assuming temperature independent 
thermal conductivities of the pure Ge and Si materials, 
justified here by a narrow AT, the heat flow in one or the 
other direction depends only on the thermal boundary 
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TABLE II. Optimized parameters Ri^aa and A/ according to 
eq. for thermal transport from Si to Ge and vice versa for 
Si/Ge interfaces with different bonndary thickness dj. 





Si Ge 



di [nm] 

Rl,oo [ 

m^K] 
GW J 

A/ 

[nm] 

0 

3.76 

± 

0.17 

167.2 

± 

13.6 

0.5 

3.66 

± 

0.12 

78.1 

± 

6.4 

1 

4.38 

± 

0.17 

58.9 

± 

7.5 

2 

7.18 

± 

0.21 

26.6 

± 

4.5 




Ge 

Si 



dj [nm] 

Ri,<x [ 

m^K] 
GW i 

A/ 

[nm] 

0 

5.76 

± 

0.28 

115.6 

± 

11.8 

0.5 

4.66 

± 

0.12 

76.3 

± 

5.0 

1 

5.70 

± 

0.08 

53.7 

± 

2.4 

2 

8.34 

± 

0.25 

27.8 

± 

4.7 


resistance (eq. [^. Following this, we define the rectifi¬ 
cation factor fuect of the TBR as the difference between 
the thermal resistance of backward (Ge to Si, Rb) and 
forward (Si to Ge, Rf) heat flow divided by the one of 
forward heat flow (eq. [Io| ). 

fuect = (10) 

Thermal rectification fuect as calculated here includes 
only the contribution of the TBR to the heat flow but 
neglects the effect of temperature dependent heat trans¬ 
port in the crystalline segments. The rectification of the 
TBR remarkably drops from 0.53 to 0.27 when the sharp 
boundary between Si and Ge is smoothed to a bound¬ 
ary layer with finite thickness of 0.5 nm (Fig. [^. Fur¬ 
ther increase of the boundary thickness only had little 
effect on the thermal rectification. This is in agreement 
with what has been shown previously for thermal recti¬ 
fication of a two-segment device. Rectification in such a 
device is increased when the asymmetry and mismatch 
of the different materials in the two segments is more 
pronounced. The introduction of a boundary layer with 
finite thickness consisting of a SijjGei-x-alloy reduces the 
asymmetry and the mismatch between the materials that 
are directly connected, thus resulting in a decrease of the 
rectification factor. 


IV. CONCLUSIONS 

Approach-to-equilibrium molecular dynamics 
(AEMD) simulations have been applied to deter¬ 
mine the thermal boundary resistance (TBR) of Si/Ge 
interfaces. The overall thermal conductivity has been 
described by a connection of series of several resistances, 
including the resistance of the interface. Its dependence 



dj [nm] 

FIG. 6. Rectification fuect of the bulk thermal boundary 
resistance of Si/Ge interfaces as a function of the interfacial 
thickness dj. It has been calculated according to eq. 

on the sample length has been described by a reciprocal 
quadratic behavior to account for non-linear effects. 
With this assumption the length dependent TBR could 
be expressed by a reciprocal linear function which con¬ 
verges to a bulk TBR for infinitely large simulation cells. 
The bulk TBR for thermal transport from Si to Ge at 
300 K resulted in 3.76T0“® m^K/W in agreement with 
previous calculations using non-equilibrium molecular 
dynamics simulations. 

The effect of the interface morphology on the bulk 
TBR has been evaluated focusing on the variation of 
the interfacial thickness. The interface, consisting of a 
Si 2 ,Gei_a; alloy, has been increased from a sharp interface 
to one with a thickness of 2 nm. TBR is found to slightly 
decrease when switching from a sharp interface to an in¬ 
terface with finite thickness. This can be explained by 
reduced mismatch and asymmetry between the pure Si 
and Ge and the interface consisting of a SiGe alloy. Fur¬ 
ther increase of the interface thickness, however, leads 
to enhanced phonon scattering and results again in an 
increase of the TBR. 

Furthermore, the effect of heat flux inversion has been 
investigated simulating thermal transport from Ge to Si. 
In this case, thermal boundary resistance has been found 
to be higher than for thermal transport from Si to Ge, 
independent of the thermal boundary thickness. From 
these results, the rectification of the Si/Ge interface has 
been determined. It is most pronounced when the inter¬ 
face is infinitely sharp (0.53) and decreases significantly 
when the interface is smoothed over a finite thickness. 

These results give insight into the thermal transport 
properties of Si/Ge interfaces, indicating a notable influ¬ 
ence of the composition and morphology of the interface 
on the thermal boundary resistance. At a certain bound¬ 
ary thickness phonon scattering dominates over the ef¬ 
fects of reduced mass mismatch at the interface, leading 
to an increase of the thermal boundary resistance which 
even exceeds the value at a sharp interface. It is thus 
suggested that reduction of thermal conductivity is more 
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effective at extended interfaces with a certain thickness. 


have been conducted on the HPC resources of CINECA 
under the project ISCRA_THETRASL 
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